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ABSTRACT 

Two new methods, suitable for either hand or digital computer 
calculations, are developed for the determination of all the trees of a 
nonoriented linear graph. Computer programs based on these methods are 
written, They greatly reduce the computation time while minimizing 
computer storage requirements. 

In addition, a program for the determination of network functions 
is also developed by using topological relationships. This program is 


sufficiently generalized to permit its use for either a one port or a 


two port passive network without mutual inductances. 
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i Introduction. 

Topological formulas for analyzing electrical networks have become 
popular in recent years. The main reasons for the current interest in 
this field are that they are relatively simple to manipulate when com- 
pared with the conventional method and that topological formulas read- 
ily implemented by means of computer programming. The topological 
formulas for analyzing an electrical network and for computing network 
functions have been studied by a number of investigators [1]-[8]. In 
section 2, topological formulas for evaluating various network functions 
of a passive network without mutual inductances are briefly discussed. 

‘Trees of a network play an important role in topological network 
analysis. How all the trees of a network can be found effectively and 
economically is the main task with which this thesis is concerned. In 
this paper two new tree-finding methods, which are named the 'dichotomy 
method' and the 'T-triangle method', have been developed. The main 
advantages of these two methods are that there is no duplication, com- 
putation time is reduced, and they are suitable for both hand calcu- 
lation and computer programming, The algorithm and computer program 
for each method are included in section 3. 

A computer program for the determination of network functions is 
presented in section 4. The topological formulas, modified by Chan and 
Chan [6] and [8], which permit the determination of network functions 
by means of only a single tree finding subroutine, were used in this 
program. This program can handle a rather complex network due to the 
smaller memory requirement achieved by the 'T-triangle' tree finding 


method. ; lf 














Dye Topological Formulas. 
2.1 General Description. 

The basic laws used to solve for the currents and voltages in an 
electrical network are Kirchhoff's current law and Kirchhoff's voltage 
law. By using these laws one can write a set of loop (mesh) equations 
or a set of node equations from a network. The ordinary method of 
solving a set of equations to obtain a desired network function from 
the ratio of a determinant and its cofactor is cumbersome. Since the 
information in the determinant contains a great deal of redundancy 
(i.e., in the expansion of the determinant many — occur more than 
once, but with opposite signs), its evaluation will be lengthy due to 
cancellations. As we see from the properties of an incidence matrix 
and circuit matrix, [1], [2], the determinant and cofactor of a system 
equations can be evaluated directly from the trees, cotrees and K- 
trees of the network. The topological analysis of electrical networks, 
therefore, provides a shortcut in evaluating network determinants and 
cofactors because the usual cancellations inherent in the evaluation 
of the network determinants are avoided. Moreover, it displays in a 
very intuitive manner, the causal relationships between the several 
variables under study, 

In network synthesis, trees of a network will give the important 
information that significantly affects the network functions, Conse- 
quently, trees of a network play an important role in topological 
analysis and synthesis. 

2.2 Definition of Terms and Symbols. 
Edge or element - a line segment together with its end points. 


Vertex or node - an end point of an edge. 
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a — a 


Graph - a finite collection of edges such that no two edges have a 
point in common that is not a vertex. This is also called a 
linear graph. 

Subgraph ~- a subset of edges of the graph. 

Connected graph - a graph G is connected if there exists a path between 
any two vertexes of the graph. 

Circuit or loop - a subgraph in which every vertex has two edges 
incident to it. 

Tree - a connected subgraph of a connected graph containing all the 
vertices of the graph but containing no circuit. 

Cotree - a set of chords, which includes all the edges of the 
complement of a tree in a graph G. 

2-tree - a subgraph of a tree of G formed by removing any one edge 
from the tree. In other words, it is a pair of unconnected 
circuitless subgraphs which together include all the vertices of 
the graph, 

3-tree - three unconnected circuitless subgraphs obtained by removing 
any one edge from a 2-tree of the graph. 

Incidence matrix, A, - a matrix with rows corresponding to vertices and 


columns to edges of a graph, and the entry 


a = 1 if the edge j is incident at vertex i 
4 = Oif the edge j is not incident at vertex i 
Circuit matrix, B. - amatrix with one row for each circuit in the 


graph, and the entry 


Bes 
tJ 


1 if the edge j is in circuit i 


aes O if the edge j is not in circuit i. 


il 








A = incidence matrix with any one dependent row deleted 





B 


circuit matrix with the dependent rows deleted. 

2.3 Concept of Topology in Electrical Networks 

Kirchhoff's Laws 

The fundamental equations of electrical network theory are: 

(a) Kirchhoff's current law - the sum of the instantaneous 
currents out of any node are zero; in matrix notation 
it is 

AI (T) = 0 | (2-3-1) 
where A is the incidence matrix of a network as 
defined in section 2.2 and T(t) is a column matrix 
of branch current. 

(b) Kirchhoff's voltage law - the sum of instantaneous 
branch voltages around any loop is zero; in matrix 
notation it is 

B VT) = 0 (2-3-2) 
where B is the circuit matrix of a network as defined 
in section 2.2 and V(t) is a column matrix of branch 
voltage, i.e., t 
Volt) = Le Fe Tet) + Ry Lett) + Ce [Tut Wojtl, © 


and where 


Le = element inductance matrix 

R, = element resistance matrix 

Cc, = element capacitance matrix it 
E(t) = column matrix of branch voltage sources 


The Laplace transformations of the fundamental equations are: 


A 1,(S) =" 0 (KCL) (2-3-3) | 








BV.(S) = 0 (KVL) (2-3-4) 
where V5 (s) = Ee (5) + Ze (Lets) — L105 + LM) 


and Z is the element impedance matrix in square form, 
= 





¥ When we discuss network functions the initial conditions 
are always set to zero, therefore 
= 2-3- 
VCS) E,(S) + Z,(S) 1(S) ( 5) 
' 
and 1,(S) = B 1S) (2-3-6) 
1 
VCS) = A V6) (2-3-7) 
where A and B are the transport of A and B respectively; 
Ih and We both are column matrices. The subscripts m and 
" n imply mesh (loop) and vertex (node) respectively. 
Substituting (2-3-5) into (2-3-4) we have 
BLE+ Ziwlsl= o (2-3-8) 
Substituting (2-3-6) into (2-3-8) we have 
oy ay Ff 
(SLA ells) = 1B Eas) (2-3-9) 
Referring to equations (2-3-5) once more, solving for Is) and 
using z.'(s) = ¥,(S)5 we obtain 
Te(s)= YO)[ Vets) — £, (5)] (2-3-10) 
where Yo is the element admittance matrix in square form. 
q Then substitute (2-3-10) and (2-3-7) into (2-3-3) which becomes 


AX OA VaS=A KOE) | (2-35-11) 


Equation (2-3-9) can be written in the form: 


z,(8) T(s) = El (s) (2-3-12) 





where Z_(s) = 8 Z(s) B ,» the mesh impedance matrix and Es 

is the column matrix of voltage sources around a loop. The 

polarities of the voltage sources are chosen to be positive 

in the direction of the loop current. 

Equation (2-3-11}) can be written in the form 
¥ 6S) V6) = I (5) (2-3-13) 

Where ¥, 6S) =A YA » the node admittance matrix and 1a repre- 
sents the node current sources in column matrix form, 

Equation (2-3-12}) is known as the loop (mesh) equations in matrix 
form and equation (2-3-13) is known as the node equation in matrix 
form. 

The utilization of the properties of incidence matrices and 
circuit matrices and the development of topological formulas to eval- " 
uate the mesh impedance and nade admittance determinants and their 
cofactors by using tree, cotree and 2-tree concepts have been dis- 
cussed by Mayeda [1], Weinberg [2], Chan [3], Reza [4], and Seshu and 


Reed [5]. Therefore, the formulas are merely outlined here. 


(1) A=| 7, (s}| = = Tree admittance products = V(Y) 
” (2-3-14) 


i.e., the sum of all tree admittance products of a network, 


N. The notation 4 is the node determinant. 


(2) Az = | A, OA. a >. Ve: = admittance 


products = Or ( Y) 





(2-3-15) 


i.e., the sum of all the 2-tree admittance products of a 
network, N, with each 2-tree formed in two subgraphs; one 


subgraph containing node i and the other subgraph con- 


taining the reference node r. This is identical to the sum 





of all tree admittance products of the graph obtained by 

identifying the go node with the reference node. This 

graph is designated N, e+ The notation As; is the (i,i) 
> 


cofactor of A . 
(3) 5 | Zn(5)| = => cotree impedance products (2-3-16) 


i.e., sum of all the chord set (complement edges of a 
tree) impedance products of a network, N. The notation A 


is the mesh determinant. 


(4) Bi: om: | (3_; ZS) 3. = = cotree impedance products 
(2-3-17) 


of N) = CLM} 


where Ny is the network obtained by deleting element y 
which is an element in loop i but not in any other loop. 
The notation As is the (i,i} cofactor of A. 
2.4 One-port Networks 
One port network (1-p network) is a network with one pair of nodes 


designated as “input terminal nodes," i.e., the driving point terminals 


A one port network is shown in Fig. 2-1, 





| 
ee | 
| 
Faq. 2-3 
i1-P network i-pP network , \—p metmork 
Rx Cited by current Rxcited by Voltage 
Aourcg, dAource 
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1 
The network can be excited across input terminals 1-1 by either 
a current source, Fig. 2-2, or by a voltage source, Fig. 2-3. 


Driving Point Admittance 


¥ 
The driving point admittance is defined as 


Tap (Ss) iL 
ie Ge ap = ‘ f. 


Vap‘s) initial condition =o V, v4 


Fig. 2-4 


1 - p network with node 1 


as the reference node. 





Starting from the node equation for a l-p network with v-1l 
independent nodes, from (2-3-13) } 
¥, 63) V6) = I (Ss) (2-4-1) 


or in detail 





Yi TR <= ee Ty, Vel Vii 1 (2-4-2) \- 

le Yow h*? Vavs Var = P 

Yeu Yuu 7. <a baer Vet’ ae) 

The solution for Via! is 
A 1 

Vii = i I, (2-4-3) | 
then we get | 
a ' 
Yap (3) — we (2-4-4) 
ta? : 


Since A and M,y are the determinant and its cofactor of 
the node admittance matrix of the network, N, therefore topological 


formulas (2-3-14) and (2-3-15) can be applied here, and the desired 


topological formula for the driving point admittance is 





ee , 
dp Wy) (2-4-5) 


. ). tree admittance products of network ey wv | 


or Y4,(S) mi ——— aye eo} 





tree admittance products of network ‘e Nay 


Example 2.1. A typical network is shown in Fig. 2-5. 





Fg 2-6 
Tha graph of natok Af, 





Table 1. The given values of parameters 





Trees of graph N; 


Lh ae a4 , 


¥n%=6S 2% = 





4 % Ye +. ae a a 


x uot 
ro 
ra 

w 

N 


¥ Bes WeGes YW =e Aa 
? i Se 
. VY) = 7) + 5 + 12 + 66 
s 
To get Wy 1! (Y), we short nodes 1 - 1' then we get Ny 1! 5 


shown in Fig. 2-7 


Fig. 2-7. The graph of Moa! : | 





Trees of graph N ; are 
Lh 


% ’ 
x 
4 rr £ 
%3 Ve Yes | 
yo 
% 
% “2 


yyezs wees 
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Wye t+ Brio + i244 bd (2-4-7) 
Therefore 
; c p+g tt? +64 — 
dp tr>~A A wvtorwa ree 
i +2 a+ ba! 
‘ap ~ i or + ao A+ 6a? a 


Driving Point Impedance 


The driving point impedance of a I-p network N, see Fig. 2-4, 
is defined as 
ahs initi 0 (2-4-10 
24,{8) = Ts) all initial values = (2-4-10) 
From equation (2-3-12), the loop equation of Il1-p network N 


with £ independent loops is 


Zin ay, 2? Cee £ie I, Ej 
fa Man oo Mot i ak | (2-4-1) 
O 
Lay hee aS ik - hy Ty 
where ae = the impedance in S-domain common to loop j and k 
t, = loop current around loop j 
Ej = net voltage rise around loop j 
Solve for I 
A 
= an E, (2-4-12) 
so that 
A 
Z = (2-413) 
ey 


Zn (mesh impedance matrix) in equation (2-3-12) is not just for 
the network N, but in addition, ineludes one loop which consists of the 


driving source, i.e., the loop established by the first row in equation 
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(2-3-12). Z for the network alone is the submatrix of Zn obtained 
™m 

by deleting the first row and first column, Therefore, with a little 

modification to the topological formulas (2-3-16) and (2-3-17), [1], 


we have 


A, = CLV @] (2-4-14) 


A ime C [ Ww, i" ( z)) (2-4-15) 


Hence the topological formula for the driving point impedance is 


C(W,1(2)] 
se ee yA (2-4-16) 
Lap C[lv @] 


To avoid the cumbersome manipulation on cotrees, we can compute 


Z from Y 





dp dp 
| | 
ver Peigeieed ee (2-417) 
M A, 3 Tap 
By the topological formula of Yap » equation (2-4-18), we have the 
more useful topological formula for His » thus 
71? 


Zan = VM (2-4-18 ) 


Example 2.2. For the purpose of comparing results and in order to 
provide the same input for a computer program, let's consider the same 


l-p network as in example 2.1. Also the given values of parameters are 


the same as in example 2.1. 





Graph of network N: z 


Graph of network N, |, Zz 
1a" Ae 4 


4 


! 
[A elena 





é 








! 
roa! zante 
i i 


Te —” ~ 





en ne RE 
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= " { 
A ee ele ya HAL 
12 12 12 
ais = <2 ae (2-4-19) 
| + 4A + 1047 412.4% C20 (24-20) / 
Z CC orl j 
ae | + 74 + 1247 464 
» 
If we use equation (2-4-17) in conjunction with the value of 
i obtained in example 2.1, then we have 


; | +4a 4j0ed' 4+ j2a' dca? 
ee ee eee (2-4,-51) 
oP Yap 1+ 74 +120°4 643 





Notice this result agrees w th the previous one. 
‘2.5 Two-port Networks 
A network with one pair of nodes (1,1') designated as the input 
terminals and one other pair of nodes (2,2') designated as the output 
terminals is called a two port network (2-P). It is shown in Fig. J 


2-8. 





Va and (2,2') desig- 
| nated as the input 
2° and output terminals 
respectively. 








i, . 1, Fig. 2-8. A two-port network, 
I K 4 t a with nodes (1,1') 


There are four variables to be related in a two-port network, i.e. 
I> the input current; Vi> the input voltage; I5> the output current; 
and Vos the output voltage, The topolotical formulas associated to 
the ratios of the different variables of a 2-p network will be discussed 
in the following sections. 


Open Circuit (0.C.) Impedance Function ks 


From Fig. 2-8, we can write the loop equations as 





Vi(A) Zu (4) Z,,(A)|} | Ta) 
= (2-5-1) 
V2 (A) Zz t (4) Xu (a) L A) 


The parameters 24> Zi9° Zo and Zoo are called the open 
circuit impedance functions of the passive two-port network. They 


are denoted as Zoe , the qpen circuit impedance matrix, thus 

Zi, Liz 2 

Zoe = (2-5-2) 
211 £32 


From Eq. (2-5-1), clearly we can evaluate these functions in the 


following manner: 











Vv. V, 
2 al tie. 
11 LW I,=0 Bik I,=0o 
(open 2-27) Copem 1-17) 
(2-5-3) 
V; V: 
Zo, =—— L.=—= 
2 iy IL=0 ie 1,=90 
(open z-2") Copenh 1-1’) 
Notice for a reciprocal network we have 
Z19 = Zo4 (2-5-4) 


Let node 1' be the reference node in'Fig. 2-8, and assume there 


are n independent nodes in the 2-p network, then the node equations 


are 
YY + %aVet cree You Va = Ty 
ead, =i You We* a es ee Yen Vn = io 


: : ; ‘ (2-5-5) 
Ter ¥, + Yee Ve Ee Ye YonVn = In 


Solving the above equation, [3] [5] we get 
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Bi Ariz — 2072’ 


| 
ence (2=5-6) 
A | Ay 7A;2’ Baz + Aa’ —Az2’/- Aa’ 


Lee seal 


where A and Aa, are the determinant and its cofactor of the node- 
admittance matrix of Eq. (2-5-5). 


From the topological formulas stated in Section 2.3, we have 


Aw om lah) ait = WY 2,77 — Wi 2/1” 
.. (Wra,1-2 +W,22' = (Wrz j2 + Wi 227 ’) 
= Wye 2! —Wee’ /e (2-5-7) 


and 


Baz + Azz! —Azrz-Axz = Ws, 1’ + W2! i Wh 2! 1 —We'e, / 

= W,, 12? Waz! 1/4. Weezy) t Wez ‘ee Wa" re Wart” 

= Wy 112? + Wet pe 

= We, 2’ (2-5-8) 
Therefore, the topological formulas aor open circuit impedance 


functions are 


1 oh Wii Wrz, 14° —Whasi’e 
Ve hee | 
oc 
vy) Wha, 72! —W 2/2 Wa, 2? (2-5-9) 


However, the more useful topological formulas of Zoe for 


digital computation are developed by Chan and Chan, [6][8], that is 


| Wh (Whe +0," -W Wa) 


Loe® 
c (Y) 5 (Wa! +Wa, 17=Wy,2-Wir2") Wie” . pene 


7 





a 


en fees 


Short Circuit (S.C.) Admittance Functions 


If we solve Eq. (2-5-1) for I in terms of Zoe? V 








1 and Vy 
; then we have 
i, oe nO “La (2-85-11) 
: =[Ze] | l= > 
iy Va ae ~Lar Bik 
We define the short circuit admittance matrix Yoc as 
a i =F Gr Na, (2-5-12) 
sc = : 
e Ya nz 
Then the entries of Yoo » which we called the short circuit 
admittance functions of the passive 2-P, can be measured in the 
following manner: 
Y= =| at ie 
ul V; V2=0 1 VV. V, =o 
; ( shovt 2-27) Csheat \-) 
(2-5-13) 
et 22 ‘= 2 
ae ag Va=e oe YW=0 
! ' one 
(shetk 2-2) (short 1-1’) 
Also notice for a reciprocal network we have 
Yio = Yor (2-5-14) 


The expressions for the short circuit admittance functions in 


terms of the cofactors of the node-admittance matrix are given (39[ 2} 
as 


| AzatAz/2'-LAY 





i 4,2! -Orz 
2-5-15) 
Sc ( 
A, )22 tA, 22 “LA y2ah Are? — Ayre 4&,, 


where As see is the determinant obtained from 4A by deleting rows 
i and k andcolumns j and 4, 
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The denominator on the right hand side of Eq. (2-5-15) can be 


expressed in topological formula as [3][1] 


Ay, 12 + By ai)\= 2A, a7 = 2 UC) 


=U, Boi. a Up 2,12’ + Ui2,21)'+ U2 2 (2-5-16) 


where U is the sum of 3-tree admittance products, The 


2g ds" 
3-tree designated by the indices of WU results from the subgraphs with 
vertices 1,2' in one part, vertex 2 in another part and vertex 1' in 
the third part. Similar interpretations can be applied to the other 
terms in (2~5-16). 

The topological expressions for those terms in the matrix of Eq. 


(2-5-15) are shown in the previous section. Consequently, the topol- 


ogical formula for short circuit admittance functions are 


| Wax Wea i’ — Wea, 12! ? 
%. =— (2-5-17) 
=U W, 27,742 — W212’ W, 1? 
However the more useful topological formulas for gc for digital 
computation are given by Chan and Chan [6] [8], that is 
| Wa,2" £ (w+ u2t ~W.2'- Wa, ¥) 
i art (2-5-18) 


= U + CW). + Wy 2—W) 2’—We, 1) W,, ‘* 


It has heen shown by Chan [7] that 
e, U = S Tree admittance product of the network which is 
obtained by identifying nodes 1-1' and nodes 2-2'. (2-5-19) 


il@hs 


=U-= WY) of the network a # (2-5-20) 
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Voltage Ratio Transfer Functions 


The voltage ratio transfer function, hig » of a 2-P Network 


(Fig. 2-9) is defined as the ratio of the output terminal voltage to 


the input terminal voltage when I, = 0. 














Fig. 2-9 
Thus 
NM, = cs (2-5-21) 
72 V, %,=06 ( open 2-2’) 
Similarly 
V, 
M=— (2-5~22) 
al Vz |r=0° (open 1-1") 
From Eq. (2-5-1) and Eq. (2-5-10) we obtain the following 
equations: 
Vv, | : Wi, 1? + (wi, 2 Wh’ -W,2-Wyz J, 
= ——— ’ (2-5-23) 
v3 v(y) = AW, at We, -Wh2-Wy; 2’) Wz 2! I, 
when terminals 2-2' are open, I, = 0. Then, 
W) ye 
V=—et I, (2-5-24) 
ve = 
= (Wi, 2’ +W —Wi,2—Wiz2') . 
\, = > ~ 1, (2-525) 





VW) 





Therefore, 


Va _ o (Wig + Way —Wi2 -Wy 2’) 





Mia = 














= 2-5-26 
V, 1,=0 Wi ( ) : 
Similarly vs 
\ 
V, = Wi, 2° +A) enw, —Wy, 2 
Ax = 1s a ( 4,1 L2 ' 2’) (2-5-27) 
V2 lee Wa, 2’ 


Current Ratio Transfer Functions. 
The current ratio transfer function, Qin» of a 2-P Network 


(See Fig. 2-9) is defined as the ratio of the output to the input 


currents when the output terminals (2-2') are shorted, i.e., Vo = 0, 
Thus 
12 
a= (2-5-28) 
42 1, lva=° C shert 2-2') 
Similarly 
A, = a (2-5-29) 
ie JT. Ivyezo chert i-v) 


From Eqs. (2-5-11) and (2-5-18) we have the following equations: 


i, Ya ie v, 
i, ¥ Yas Yor V: 
| Wa, 2’ 4 (Wat Wiz! -Wi2'-Wd VY 


* 2.3 4(W, 2+ Wy2’ Wi, 2’-Wa,v) Wo V.K2-5-30) 
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By similar manipulations to those done for the voltage transfer E 


function, we obtain 





Ms 4(Wia +Wi2! — Wy, 27 —Wa2,v) (2-5-31) 
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Given a 2-p Network shown in Fig. 2-10 
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(b) The graph corresponding 





(a) 

Fig. 2-10, (a) A given 2-P Network, 
to the Network shown in (a). 

The values of the parameters in the Network are shown in Table 2. 


Table 2 Parameters of Example 2-3 
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The network functions of the 2-P Network given above are obtained 


by the following three steps: 





Step 1. List of all the trees, 2-trees, 3-trees and their admittance 


product. 
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Step 2 


From the results of Step 1, we have 


VY) = a4 E4444 A° 

WpyWs= 4244040 
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Step 3. By the topological formulas of (2-5-10), (2-5-18), (2-5-26), 
(2-5-27), (2-5-31) and (2-5-32), we obtain the following results. 


(1) Open circuit impedance functions 
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(2) Short circuit admittance functions 
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3:5 Tree Finding Method 
Most topological formulas for network analysis depend on the 

determination of all the trees of the network graph, The number of 

trees of a complex network may be greater than a million. Therefore, 

the computation time and memory space are the primary factors we have 

to consider when a computer is used for tree finding. Many tree-finding 

methods have been given by various authors [9] - [15]. But one or 


more disadvantages exists in each of these methods, They are listed 


in Table 3, 














TABLE 3 

Tree-finding method given by Disadvantages 
G. J. Minty i long computation time, since 
(2T-1) graphs are needed for 
j T trees. 

G. W. Zobrist and simplify the tree finding method 

G. V. Lago for a complete graph only by one 

: step. : 

| i 
i W. Mayeda and large memory space required, also 
j S. Seshu tedious for hand calculation, 
i 
M. Piekassky long computation time and Large 
; memory are required due to 
; duplication, 
Wai-Kai Chen Same as above. 
: , } 

B. R. Myers and ' Same as above. 
; L. V. Auth 
L. M. Maxwell and Same as above. ‘ 
i J. M. Cline 





Two different tree-finding methods are given in this paper: 
Dichotomy method and T-triangle cebhod., Both of these methods do not 
result in duplications, hence, these two methods achieve both a reduced 
computation time and a smaller memory requirements. However, the T- 
triangle method is much better for a complicated network, 

3.1 Dichotomy Method 

Algorithm 

1. List the sets of edges incident to each of the vertices of the 
given graph and arrange them in any arbitrary sequences in a row 
(see examples) and call them "vertex sets." 

2. Comparison Procedure 
(1) Take the first vertex as an "operating vertex." Compare its 


vertex set with one of any other vertex which has at least 
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one edge in common with that of the operating vertex. This 
other vertex is designated as ‘comparing vertex’, List 
(store) the common edges found by this step as a “common set." 
(2) Create two new "vertex sets" in the following manner: 
(a) Eliminate the common edges from both the operating vertex S 
set and the comparing vertex set. Place them with the 
remaining ordered vertex sets to the right of the same 
row. 
(b) Combine the operating vertex set and the comparing vertex 
set with their common edges eliminated to form a new 
"operating vertex" set, then place it together with the 
other original vertex sets below in the same column, 
3. Repeat 2 successively for the new vertex sets through the same « 
column until an empty set is resulted, 
4, Return to the first non-zero vertex set, which was stored in the 
preceding rows, and repeat 2 and 3. Repeat this successively 
until an empty vertex set is found in the first row. : 
5. List all the trees by the following steps: 
(1) Form an ordered group from all "common sets" in the first 
column. Form additional ordered groups by taking the common 
sets generated in each succeeding column and replacing them 
for the common sets in their respective rows in the preceding 
column until the common sets of the final column have been 
used, 


(2) List all the trees from each group of common sets by taking 


one element at a time from each of the common sets. 





Example 3.1 


Fig. 3-1 The given graph with each 
edge designated by a 
number, 





(I) List the vertex sets 
[12] (134) (2356) (457) (67) 
Operating set. 
(II) Comparison 
[Legh Seas Seats HET) em @) 
J 


[2 34}(23 56) (45 VET) 
* * xe 





—w (4}(s6) (4576 7) 
* 


23 4 
{4 5é\(45 T)(6T) wm [6] (77) LST) 5 6) (67) we [7] (6) (67) 
x * ud « « * & a * 
RS z, = 7 
6 
lea" NEIGH C67) (¢ 7) [6] 6) 
67 T 67 G 
@) [.2\(34)(2356)457)(67) ——________+ [0] --- - step 
« 3 
z 
[354] (3 A)C45 7) (67) ———» U6 ]ea) 4s 797) [6] )(47) (67) 
¥ % Pa 
3 $ 6 
[5 64) (457)(67) ——m[6}(7)(67)  [647)(4)C6 7) C7] (4) ¢47) 
et oo RE % * * * om * 
45 é * 7 
[67] (67) (747) Lé 7] (67) Ce} c4) 
ee ME = om xt 8 ® * % 


67 e a 4 








(III) Listing all the trees 


Group of Common Sets TREES 


12 AG 1256 13426 136 
(1)(23)4 5967) |, 2,47  sags7 1347 meet 


(V)C23)4e>(7) 1267 13987 


Gdyc 470 S (eT [vars /4a57 


Cha O29 7 NG). eZ 


(2)€3)(45)G@7) | 2246 232467 2366 2367 


(2) € 37H 6)C7) | 2367 


(2) 6 S)(4)(67) | 26°46 2547 


(2) (6) (7) (4) | 2474 


Example 3.2. 


Fig. 3-2 The given graph with 
each edge designated 
by a number, 





(1) List the vertex sets 


[123] (145) (36) (2456) 
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(II) Comparison 


[123] (1453 6)( 2456) ———_____—» [2 3] (45 )(36)(24$6) +@) 
% *> * * 


" | 3 
(2 345)(3 6)( 2446) [245] (6X2456) (2€](45) (2456) 
* * w tm TER « % ik m 
. 3 245 Zé 
{2456](245 6) £éJ(6) [45] (45) 
mK E* RK x «© *® ae x a <a 
2458 6 45 


@ pad er LOO C PSE &) —r Lol stop 


= 
[456](45)(6) 
* % » *% 
45 
7 
[e}(é) 
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(III) List all the trees 


Group of Common Sets 





Total number of trees = 13 
Discussion 
As illustrated by these two examples, we have the following con- 
clusions to this tree-finding method. 
(1) All the combinations which listed from each group of the common 
sets by procedure 5, are trees, 


(2) No duplication of trees occurs, 
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(3) All the trees of a given graph are listed. 


The proof of these three statements are given below: 


(a) 


(b) 


The first common set are the edges common to the operating 
vertex and the comparing vertex, in other words, the corre- 
sponding edges in the first common set are incident to the 
two vertices of the graph. The elements in the second 
common set correspond to the edges, which are common to a | 
new vertex and one of the prior two vertices (the vertices 
had been compared for obtaining the first common set), 
Similarly, the associated edges of the elements in every 
additional common set will connect a new additional vertex. 
Since there are V-l (V is the number of vertices of a graph) 
common sets in each group, a combination, taking one element 
trom each common set of one group, will connect V different 
vertices (all the vertices in a graph). Therefore, the 
edges collected by such combination is a tree. 
J W mw! . ° th 
Since the new ''vertex sets'' stored in the i+l column are 
generated by eliminating the common edges which have been 
F .th 
found in the i column of the same row (see procedure 2-(2)- 
: .,,th 
(a)), the elements in the common sets of the i+l column 
will exclude some elements which are in one common set of the 
.th : F 
i column. Consequently, no trees list from different group 
of common sets will be the same. Thus, there is no tree 
duplications. 


We start our procedure by taking the first vertex of the 


original vertex sets to be the operating vertex; then we 





follow successive comparison procedure to find the common sets. This 
means that we take the edge associated to the elements in the operating 
vertex (say vertex A) as the reference edges and trace through them to 
the other vertices. The possible ways to reach other vertex (say 
vertex B) ies through the edges which are common to both vertex A and 
vertex B. Other possible connections involve the edges corresponding 
to the elements, which are resulted from the ring sum operation of the 
vertex sets of A and B. This set of elements constitutes the new 
operating vertex by procedure 2-(2)-(b). The following Figures illus- 


trate this point. 


Ba 


(a) Operating vertex A (b) Common set ANB = hs 
New operating vertex 
set: ABB = [2,3,4} 


As we continue this process until one vertex set is empty, we get 
one group of common sets. Moreover, in each comparison step we also 
generate another new 'vertex set' with those elements of the common 
set (which have been found " the previous column) eliminated. By 
repeating this process repeatedly, we have all the possible groups of 
common sets. Thus, all the trees of the graph are found. 

3.2 Tree Finding Computer Program I (by dichotomy method) 


A digital computer program has been written in FORTRAN IV language 


for finding trees by the "dichotomy method." This program is included 





in Appendix I. Examples 3.1, 3.2 and 4-vertex to 8-vertex complete 


graphs are solved by the same program, General descriptions of the 


program are given below: 


(1) 


(2) 


Inputs of 
the second 
problem 


Inputs of 
the first 
problem 


Program Limitations. 


(a) Nonoriented linear graph. 
(b) Maximum number of vertices of a given graph is 15. 


(c) Maximum number of edges in a given graph is 60. 


Input Data. 


























7 Second row of an incidence matrix 
First row of an incidence matrix 
(No column space) 





Incidence 
matrix 


(KAM) 






MB 
ft 2! 134 


Fig. 3-3 Input Data Cards 


N - Number of vertices of a given graph. 


MB 


- Number of edges of a given graph. 


LIST - Conditions for tree printing. 


LIST 
LIST 
KAM 


(3) 


= 1: print out all the trees that have been found. 


0: do not print out the trees. 


-~ the incidence matrix of a given graph. 


Outputs. 


(a) Trees, designated as 'IANS' in the program, can be printed 
out in the form of integer numbers which are labeled for 
each edge of the given graph. 


(b) Total number of trees, designated as 'KSUM' in the program. 


(c) Program time, designated as 'CL' in the program, 











(4) Block diagram of the program. 


Set program time ~— (100) 
to zero 


Block A ' 






Read input: N, MB, LIST 
and KAM 


Initialize constants 





Block B ’ 


Print out all the input data 


Block C 






1. Convert the input incidence matrix (KAM) to 3- 
dimension matrix (KA). 
2. Set the operating vertex (KS) 









Block D a. 


















Successively call 

subroutine 'NPART' ie Pe ‘ 
for each comparison 1. Find the 'common set' | 

step. from ‘operating ver- | 


tex set' and ‘com- 
paring vertex set',. 






rt |- - 





when a group of 'common' 
1 





{ Call subroutine 'TREE' 


sets’ has been found. 









2. Generate two new 'ver 
tex sets', one for 
the same row opera- 
tion, the other for 
the next row opera- 
tion. 
























Form trees from a group 
of 'common sets', desig- 
nated as MTREE in the 
program, by taking one 
element from each row of 
"MTREE' at a time. 





















Determine the control 
conditions for the 
requirement of next 
Operation. 


















2. Print the tree out as 
soon as it has been form- 
ed when 'LIST' # 0 














1. Print out the total number 


of trees, 









2. Print the program time. 










If more than one 
problem 






SS ee) Bit ne a 





3.3  T-triangle Method. 

This tree finding method is based on a T-triangle which is 
obtained from a given graph, We generate succeeding new T-triangles 
by the procedure, which will be stated in the algorithm. Then we list 
all the trees from each of the T-triangles. This method achieves both 
a reduced computation time and a smaller memory requirement for com- 
puter implementation. It is also suited to hand computation because 
no duplications occur and no new graphs are required. 

T-triangle. A T-triangle of a graph is a triangular array, of which 
the element my > in the ‘thas row and en column of the T-triangle, 


is given by 


a : i= 1,2,-°7°.¥" 
ys ={T AG] 5S) \ 


where I is the set of all the edges incident to vertex K and V is the 


number of vertices in the graph. Fig. 3-4 illustrates how the T- 


triangle of a graph is obtained. 
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Fig. 3-4. Example illustrating T-triangles (a) a given graph G, 
(b) T-triangular diagram corresponding to G, (c) general form of a 
T-triangle corresponding to a 4-vertex graph, and (d) the T-triangle 
simplified from (b). 

If we expand the vertices of a graph to lines, called vertex 
lines, and reduce all the edges between each pair of vertices to small 
circles and connect them to their corresponding vertex lines in a 
triangular form, we will obtain a diagram as shown in Fig. 3-4(b). > 
Actually, Fig. 3-4(d) is a simplified form of Fig. 3-4(b). We can 
consider each circle, which includes all the parallel edges between a 
pair of vertices, as a single switch. When we say that the switch is 
non-empty, we imply that one or more elements (edges) exists in a 
circle. When we say that the switch is empty, we imply that no element 
exists in the circle. With this concept we can examine some properties 
of T-triangles corresponding to a connected graph with V vertices. 
Property 3.1. 

Both the number of rows and the number of columns of a T-triangle 
are (V-1L). 

Property 3.2. 


The number of vertex lines of a T-triangle is V and the total 


number of entries (circles) in a T-triangle is V(V-1)/2. 





+i 


Property 3.3 
The switch corresponding to a) of a T-triangle has one end 

connected to vertex line (i+1) and the other end connected to a vertex 

line, which appears prior to the G+ig™ row in the triangular diagram 

(or the T-triangle). 

Property 3.4 
No two switches (sets of edges incident to two different pairs of 

vertices) of a T-triangle connect to the same pair of vertex lines. 

Proof: From the definition of tay (an entry of a T-triangle), there 

is at least one subscript, either i or j, which is different from the 

corresponding subscript of any other entry. 

Property 3.5 

(a) All the switches corresponding to the entries of ‘the first column 
of a T-triangle are connected to vertex l. 

(b) All the switches corresponding to the entries of the last row of 
a T-triangle are connected to vertex V. 

(c) All the switches corresponding to the entries in the = Tow and 
the entries in the ay column of the T-triangle are connected 
to vertex (itl). 

Proof: The subscript j of each fry in the first column is 1; this 

implies that the corresponding switch has one end connected to vertex 

1. The subscript i of each “4 in the last row, the eo row, is 

V-1l. Therefore, the corresponding switch has one end connected to 

vertex V. Each of the entries of the Pe row has a subscript i. 

This implies that the corresponding switch is connected to vertex 


(i+1). Similarly, a switch corresponding to an entry in the esa 


column is connected to vertex (i+1). 











Property 3.6 


A T-triangle represents a linear graph and any linear graph can 
be converted to a T-triangle. 

The proof of this property follows from the definition of a 
T-triangle and the concepts of vertex lines and circles (switches) 
and property 3.5. 

Definition 3.1; Simple Loop 

A simple loop is a closed loop formed by two edges which are 
incident to the same pair of vertices of a graph. 

Convention: The elements inside a switch shall not include the 
vertices to which they are incident in the graph. A switch is clcsed 
to vertex k implies that the elements inside the switch are connected 
to vertex k., « 
Theorem 3.1 

When a switch is closed to its associated vertices, 

(a) No loop will result if there is only one element (edge) in that 
switch, 
(b) o simple loops will exist if there are M elements in that switch, 

where M > 2. 

The proof is sufficiently obvious to be omitted. 

Theorem 3.2 

No loop can result if only two switches in a T-triangle are closed 
to their associated vertices except those simple loops formed within 
the individual switches. 

Proof: By property 3.4 we know that two switches will connect at 


least three different vertex lines, therefore, the corresponding sub- 


graph is a subtree. 





a 
* 


Theorem 3.3 
If only one switch of each row preceding row i has been closed 
to its associated vertices, then 
(a) No loop occurs except the simple loops when any one switch of the 
Pi row is closed to its associated vertices. 
(b) Loops will result when two or more switches of the os row are 
closed to their sexdeuated vertices. 
Proof: In the first row there are two vertex lines and the inclusion 
of each succeeding row increased the number of vertex lines by one, 
Thus, there are (i+1) vertex lines when the eo row is included. From 
properties 3.3 and 3.4, the (i+l) vertex lines are only connected by 
i switches. Therefore, no loop exists, except possible for simple 
loops. If two or more switches are closed in the "he row, then (i+1) 
vertex lines are connected by at least (i+1) switches, hence one or 
more loops exist. 
Theorem 3.4 
The collection of the elements, taking one element (edge) from 
each row of a T-triangle at a time, is a complete set of branches of 
a tree. 
Proof: An element in row i represents an edge which connects vertex 
i+] and one of the vertices preceding i+l. By taking one element from 
each row, each set contains V-1 edges, which connect V vertices. Thus, 
the theorem is true. 


For example, the trees listed from the T-triangle in Fig. 3-4(d) 





hy by this theorem are shown in Table 4. 
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Table 4 Trees list from the T-triangle in Fig. 3-4(d) 






Subgraphs 






Theorem 3.5 

If all the entries of row i (for i = 2) and those directly below 
row i down to the bottom row are empty then the graph corresponding to 
the T-triangle is disconnected. 
Proof: By property 3.3 the edges which can be connected to one of the 
vertices preceding vertex i+ 1 are those elements in the rh row or 
directly below them. Since all of these elements are empty, the graph 
is disconnected. 


An example for this theorem is illustrated in Fig. 3-5. 


1 5 

a F - 

+e: 4 a ‘ 

fe) 

o.0 967N 2 ; 6 
(b) 


Fig. 3-5 Illustration for theorem 3.5 
Theorem 3.6 


If a T-triangle, T is generated from an original T-triangle, 


1? 
T , according to steps (a) through (d), 
(a) Move all the entries of the ha row of Ts directly down to the 


bottom row, and 


(b) Shift all the entries, which are located directly below the ee 


row, upward one location, and 
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(c) Shift all the entries, which are: located to the right of column 
(i+1), upward one location and toward the left by one location, 
u and 
(d) The entries in estimnetirl) are translated, successively, down- 
ward to the bottom row and to the extreme right of that row. 
(Note that all the entries prior to the a row are unchanged.) 


The graph G, corresponding to To are 2-isomorphic. 


if 


Proof: Since the procedure for the generation of the new T-triangle 
results in the rearrangement of vertex representation in the new 


T-triangle, graph G, and G, are identical with the exception of 


il 


different vertex sequences, Thus, they are 2-isomorphic. 


The illustrations of this theorem are shown in Fig. 3-6. 





Fig. 3-6 Illustration for Theorem 3.6 
Algorithm 
The procedure for tree finding by the T-triangle method is given 
below. 
' (1) Obtain an original T-triangle, Ty from the given graph. 
! (2) Take any column (call it column j) of Ty as an "operating column," 


provided that 
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(3) 


(4) 


(5) 


(6) 


(7) 


(a) the entries of that column are not all zero, or 

(b) the entries to the left of column j and below row (j-1) 
are not all zero. 

Obtain a new T-triangle, Ty from Ty by 

(a) removing the elements of the iynay ** row of To» and 

(b) shifting all the entries, which are located at the left of 
column j and below row (j-1), upward one location, and 

(c) shifting all the entries at the right hand side of column 
j upward one location and toward the left by one location. 
Simultaneously the entries of the column j are translated, 
successively, downward to the bottom row and to the 
extreme right end of the T-triangle. 

Repeat (2) and (3) for every j = 2, 3, .... V-1, where V is the 

number of vertices of the graph. 

From each of the newly obtained T-triangle (call them Ty j= 

2, 3... , V-1) from steps (3) and (4) generate succeeding new 

T-triangle (call them Tay? k= j, jtl, ... V-1) by taking an 

operating column k each time as in steps (2) and (3). 

Generate all the T-triangles successively by repeating (5) on each 

newly generated T-triangle. 

List all the trees from each T-triangle by obtaining all possible 

combinations, taking one element at a time from each row of a 

T-triangle. 


The illustrations for each step stated in the above procedure are 


given in Fig. 3-7 through Fig. 3-10, 











(a) 
Fig. 3-7 Illustration for Step (1) 





Fig. 3-8 Illustration for Step (2) 


In either (a) or (b) column 3 cannot be the operating column. 





Fig. 3-9 [Illustration for Step (3) 


53 

















Fig. 3-10 Illustration for Steps (4) - (6) 

Step (7) in the procedure is the same process previously discussed 
in theorem 3.4, 

Occasionally we may encounter a special case that all the entrees ” 
in one row of Ty or in any new T-triangle are empty. (This case 
occurs when a given graph has a large number of vertical but with many 
pairs of which are unconnected.) If such a case occurs, in the above 
procedure, we must move this empty row immediately down to the. bottom 
row as provided by theorem 3.6. Continue this procedure until no row 
is empty, then continue the other steps. This will be demonstrated in 


example 3,5. 
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Example 353 





(4,5) 6 





(a) (b) 
| Fig. 3-11 (a) A given graph, (b) Ty the T-triangle 


corresponding to the graph in (a). 


(I) Generate all T-triangles 


(3) (4) 
\ 
(4,5 
60 
ae . 





(II) List the trees of each T-triangle 


123 136 134 (= 
(1) £124 (2) 146 (3) < 135 (4), 356 
125 156 164 
126 165 
Total number of trees = 13 
Example 3.4 
2 4 4 
| 
\ / 7 Zz 3 
o 4 $5 
o oO 6 TF 
i 1 Zz 3 é 3 
(b) 


(a) 


Fig. 3-12 (a) A given graph, (b) Ty 
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(1) Generate all T-triangles 


(1) 

















(II) List the trees of each T-triangle 


1246 
(1) ¢ 1247 
1256 
1257 


2563 
(5) § 2564 
2573 
2574 


1346 1267 1475 
1347 (2) {1367 (3) {1476 
1356 
1357 


2637 (7) {2674 
(6) {3634 


Total number of trees = 21 
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| (6) 
: 


(4) 


{ 1456 


(8) { 2347 





Example 3.5 











z /# a re 
: _ ‘ eo & & a 
(a) Graph G (b): To corresponding to G 


Fig. 3-13 


(I) Generate all T-triangles 





(II) List the trees of each T-triangle 


(1) ae (2) 1254 (3) 1345 (4) 2543 
1234 
Total number of trees = 5 


ay 





Proof of the validity of the tree finding method by T-triangle 


(1) Definitions 


V 
e. 


L 


E; = 


The number of vertices of a given graph. 

Any one element in the “ row of a T-triangle. 

All the elements in the — row of a T-triangle, i.e., 
E, = { e, ie 

The set of all the trees of a given graph. 

A set at combination of V-1 elements including one 
element in row i and one element in row j of a 
T-triangle. | 


A sub set of T in which all the trees include 


elements e, and eae Also it is a subset of T'. allt 


ij 


A subset of Teze; » of which the combination of : 
V-1 elements is not a tree. That is: 


sbi = T Ut : 

ih nh ¢. 
A set of combination of V-1 elements including one 
element in row i but excluding all the elements in 
row j of a T-triangle. 
A subset of T in which all the trees must include 
the element e; but éxclude all the elements ie we : 
A subset of : -. » of which the combinations of 

i 

V-1 elements is not a tree. That is; 


oe =, 
All the possible combinations forms by taking k 


elements at a time from the set E, . Where k 


W 
= 
v 
ie) 
v 
‘ 


+e». up to the total number of elements in ES 4 
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(2) Proof: 

All the trees of a graph can be sorted out into separate groups 
according to the elements from which they are formed. If we sort T 
into two subsets such that all the trees containing element '"b" form 
one subset and all the trees which do not contain element "b" form 
another subset, then T is the sum of the trees of these two subsets. 
We shall sort the trees of the graph by recognizing the elements 
included in each row of a T-triangle. For simplicity, we assume that 
every entry, Sia , of a T-triangle has no more than one element in 
our proof, Notice that no generality is lost by this assumption, 
because when listing the trees from a T-triangle, we can consider the 
possibility of parallel edges in each entry of a T-triangle. We 


observe that 


T = Tce)t Tee, (3-3-1) 


But = as { te, | 
Thus (E,) = e, (3-3-2) 
substituting (3-3-2) into (3-3-1) we have 

T + ig) 1, (3-3-3) 
Next, Tt. can be expressed as 

Te, = Teceyt Ve~e2 - | (3-3-4) 
Notice (Es) |tuy Fan» has toa } » So thet 


? , , , 
Te (e), > Tety@ta Vete~te, + Ve tar Tar 


= (Tet, te + Te, Tay tis) * _ tat Teta-tes)* Testate (3-3-5) 
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The last term in (3-3-5) contains no tree at all according to 


7 * 
theorem 3.3, hence T, » and we can rewrite 


Eft, tar S tjten 
(3-3-5) as: : 
Te (E:) = Te, tay eet Tete ~tz + Te, (E,) (3-3-6) 
By definition we can write 7 
Tee, = Teta rts t Testis tz, . (3-3-7) 
Substituting (3-3-7) into (3-3-6), we have 
/ —_ 
= 3-3-8 
Te, (E,) Taye, aa Te, (Ez) ( ) 
Ud — eee 
es Teg) = Tete,)t Tee,» (3-3-9) 


Comparing (3-3-8) with (3-3-9), we obtain 


Tee) = Tee, (3-3-10) 
Substituting (3-3-10) into (3-3-4), we have 
Continue sorting qT, gi > we get 

=: (3-3-12 

jae = Tee) t Tee,-z, ( , 
By similar manipulations we can drive 

Te,e,c&;) = Te, 2,e, ary 
Substituting (3-3-13) into (3-3-12), we have 

Te,e, = leee,t Te,e,-£; pro) 


Substituting (3-3-14) into (3-3-11) and then into (3-3-3), we have 
T= Tage, i Te,e.-e, 7 Te -£, + Te (3-3-15) 
If we continue to sort the first term on the right hand side of 


equation (3-3-15) by the same manipulations until E. the elements 


= 5 


of the last row of the T-triangle, have béen considered, we will have 


T = Te é, wg Cy-) _ Tae, es * Cys" Eyes 5 Te, 2° -€y-3 SE es a ies : Cy4 ~Ey.3 


to: ee Teele, + Ve, -c,+ Toe e386) 





Now, let us examine each term on the right hand side of (3-3-16). 
We see that the first term gives the trees that we can list from Ty 
by step (7) as stated in the algorithm. The second term is an empty 
set (no complete tree exists) because of removal of all the elements 
of the last row (corresponding to the removal of a node from the graph; 
see property 3.5). The third term gives the trees that we can list 


ao 


-— is defined in step (5) of the algorithm.) The 


fourth term gives the trees that we can obtain from the T-triangle, 


Tt 2? which is obtained from T by removing the elements in row 
zh fo) 


(V-3) and keeping the elements in rows 1, 2, ... V-4& unchanged (in 
order that we can have one element from each af these rows). Now, with 
all the entries of row (V-3), removed, we take column j = V-2 as the 


operating column and following step (3) of the algorithm to generate 


the new T-triangle, T, , - The trees in the set ca 


2,€,°°*Cy-g-Ey—3 


are all the trees given by Fy-2° Similarly, we can sort the set 


associated with a and the results will 


2 


be the trees that we cam List from TY.» Tyg yor? Ilyas veo and 
. J 


Te, e2- vee Cy Eyig 


Te2v-2,v-1 by step (7) of the algorithm, 
By repeating the same manipulation for the other terms in 
equation (3-3-15), we will have 
The trees given by the original T-triangle, Ty > and 
all the new T-triangles, which are generated by steps 
(2) - (6), with tree listing method given in step (7) 
.in the algorithm, 


This completes the proof, 
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3.4 Tree Finding Computer Program II (by T-triangle method) 

A digital computer program has been written in FORTRAN IV 
language for finding trees by T-triangle method. This program is 
included in Appendix II. Example 3.3, 3.4, and the 4-vertex to 8- 
vertex complete graphs are solved by this program. The general de- ‘ 
scriptions of this ptagtem sire given below. 

(1) Program limitations. 

(a) non-oriented linear graph 
(b) maximum number of vertices of a given graph is 15 
(c) maximum number of edges in a given graph is 300. 


(2) Input data, 





















Inputs of df 4 ee 
the second 6 
problem, a a 
gE a news } "LINK 
Inputs of ——— nn ———— / TRIANGLE" 
the first Y First row of 'LINK TRIANGLE’ 
problem. S ee, a 
N MB LIST 
1 10| 11 20:21 30 
1 
Fig. 3-14 Input data cards 
N - Number of vertices of a given graph. 
MB - Number of edges of a given graph. 
LIST - Conditions for tree printing. 
LIST = 1: print out all the trees that have been found, s 
LIST = 0: do not print out the trees. 


The "Link triangle" of a graph is a triangular array, the element 


of which, “4 >» in the i row and ha column of the "Link triangle" 


is given by 
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be 


£.. = the number of edges between vertices i and (i+tj). 
Where i= 1, 2, ... V-1, and V is the number of vertices of the graph, 
jtisv 
an example of the "Link triangle" corresponding to Fig. 3-11(b) is 


shown in Fig. 3-15(a). 


‘a Verte x ——-} 
x 3 «O44 





nana (b) 
Fig. 3-15 (a) ‘Link Triangle' corresponding to Fig. 3-11(b) 
(b) T-triangle generated from (a) by the computer 
program. 

From the input of "LINK" the program itself will generate the 
original T-triangle, To ; and label the elements in z, by integers 
in the sequence order as shown in Fig. 3-k5(b). Notice that there 
are two edges between vertices 2 and 4 which are indicated by hoo . 
(i.e., edges #4 and #5) but edge #5 disappears in To: 

However, the computer will reeover those parallel edges which are 
temporarily omitted from the representation when the trees are listed 
from the T-triangle later. 

(3) Output. 

(a) Trees, designated as 'ITREE' in the program, can be printed 

out in the form of integers which are used for labeling the 


edges of the given graph. 


* edge #5 is embedded in #4. 
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(b) Total number of trees. 


(c) Program time. 


(4) Block diagram of the program. A 






START 


Set program time 
to zero 


Block A 


(1) Read input: N, MB, LIST 
and LINK 

(2) Convert LINK to NT (original 
T-triangle, To ) 













(3) Register the embedded edges. 


Block B 


(1) Print out all the input data 
(2) Initialize constants 


























Block C | 

oer Successively call subroutine | 
wo | "CONTL" at each row operation 
“eo | (1) Check the requirements 
5 | for the operating column > 
va | (2) Check all the entries in the ey 
: ,th tt " : S 
eS | i row for "empty or not aa aoe 
Be } (3) Recover the embedded edges ready ii 
vg a | 
ra a! for tree listing. All elements iw 
4 8 | stored in matrix 'MTREE'. e) 
= ; SS en ae ae ee ees [> 
a @.| Block D a be 
5 | Call Subroutine "TREE": ‘c 
me) (1) Form trees by taking one % 
9 @ | element from each row of a 
pea matrix 'MTREE' at a time Ft 

(2) Print the tree out as soon as yo 

7 it has been formed if 'LIST' # 0 [m4 

ee ee ee en eee wo 
Call subroutine 'NEWTRI' to generate _ : : 


a new T[-triangle 


© : 
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Block E 
he 






Print out the total number 
of trees. 
Print the program time. 










ie 





4. Computer Program for the Determination of Network Functions. 

A computer program for the determination of network functions, 
i.e., driving point admittance, driving point impedance, open circuit 
impedance, short circuit admittance, voltage ratio transfer function 
and current transfer function, has been written in FORTRAN IV 
languate and it is included as Appendix III. The T-triangle method 
of tree finding, which has been already described in section 3.3, is 
inserted in this program as a principal subprogram. The dichotomy 
method of tree finding is easier to handle a network which has 
parallel edges between any pair of vertices; but, it requires more 
memory space than the T-triangle method does. Also, the time of com- 
putation for the former will be a little longer if the network is more 
complex. The same networks given in examples 2.1 and 2.3 were solved 
using this program. 

4.1 Program limitation 
(1) Passive network without mutual inductance 
(2) The number of vertices of a given network are limited from 


Seto yt las ee 
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(3) Maximum number of edges of a given network is 300. 
(4) Maximum number of parallel edges between any pair of 
vertices is 5. 


4.2 Rules for the labelling of a given network 





(1) The labels designated for the vertices of the terminal 


ports are shown in Fig. 4.1. 





Fig. 4-1 (a) 1-P (b) 2-P (c) 2-p, when vertices 
3 and 4 are shorted together. 


The remaining vertices are labeled with the numbers shown in 
parenthesis in Fig. 4-1, where V is the number of vertices of the 
given network, 

(2) The labels of edges are designated by the integers from 

1 through E, the total number of edges. The parallel edges 
between a pair of vertices must be numbered in sequence 
order. 

4,3 Input data. 

The data card configurations required for the running of the 
program are described below. 


Data Card 


Sequence No, Columns Description and Specifications 
1 1 Blank, Dash and Plus Sign; 
2 used for printing the 
3 dividing line and signs of 


the output. 





Sequence No. Columns Description and Specifications 


4-72 Disregarded, 
| 2 1-60 Titles. 
61-72 Disregarded. 
7 3 1-50 Titles. 
51-72 Disregard. 
& 1-60 Titles. 
61-72 Disregard. 
5 1-5 Number of vertices; Max: 15 
6-10 Number of edges; Max: 300 
11-15 Ground conditions; 


1 if Fig. 4.1(c) is met. 
2 if Fig. 4.1(b) is met. 
(No constrain for 1-P) 


16-20 Port conditions; 
ae Gb 

5 2: 2-P 

21-72 Disregarded. 

6 1-10 Same integer as labeled.for the edge 
in a given network. 
11-13 Vertices connected by the 
14-16 edge as indicated in column 1-10, 


The vertex number designated in 
column 14-16 must be the larger one. 


17-21 Same integer as designated for the 
corresponding entry of the original 
T-triangle. 
22-26 Number of parallel edges between the 
r vertices as given in columns 11-16. 
27-31 Sequence number of the edge between 
4 its parallel edges. 
32-41 Magnitude of the admittance of the 


edge; in decimal point form (F10.4). 
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Sequence No. Columns Description and Specifications 


42-46 Degree in S-domain of the admittance 
' of the edge. 
47-72 Disregarded. 


The configurations of the remaining data cards of a problem are 
the same as card No. 6 and each card specifies one different edge. 


The input data cards for examples 2.1 and 2.3 are shown in Table 5. 


iE 2 3 
THtle Input data for Input data for 
(Cards 1-5 are the Example 2.1 Example 2.3 


same for any given LO kg 2 
network.) 








! ' 
1 1 B18 ‘41 51 
op —— —— 
-DRYG Y= Zil == Z12=Z21= 222 = Yil = y12=yY21= 
Y22. = V2/Vi =  Vi/v2 = -12/Tl = ~—=1/i2 = — 
vcy) W1,1'(Y) W2,2'(Y) W1,2(Y WE" (Y) WL2° te 

















Table 5 Input data for example 2.1 and 2.3 
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4.4 Qutputs. 


(1) Trees and tree admittance products of the given network 
; and modified network will be printed out for the purpose of checking. 
(2) All the network functions of 1-P or 2-P will be printed 
out in the regular form as shown in the following example: 


Z0S*#O FG.0SR HI TAOS REZE LOSERS 
DRYG Y = 
4.088##90 + 2.08% *% 1 +4.65** 241.0 Sue 


The integer which follows "**" is the power of "S,"' 


4.5 Block diagram of the program 





Block A 






T (1) Read: Signs and titles 
(2) Print out: Topological formulas which 


7 are used in this program. 





1100 
Block B 





(1) Set program time to zero. 





(2) Read input: All the information of a 
I given network, 
' (3) Form the original T-triangle from the 
| input information. ' | 
| (4) Register the embedded edges. 

(5) Print out all the input data and the 
original T-triangle, "NT." 


Block C | 


| Set constants by the ground condition, 








i 








i "IGR," and the port condition, "IPORT." 


! 
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Block D = iat. en 


(1) Set values of i and ihe 
(2) Call Subroutine "MODNET" to generage a modified 





network, J, j by identifying vertices i and 
i, > 


j of the given BEUREE Es Ne 


Block E ' 





oa (1) “Call Subroutine "TREE" to generate all the 


| 

| 

| 

trees of the network N or the modified | 

network N. . . 

a ai } 

(2) Gall subroutine "TREEPR" (a subroutine of 

"TREE') to calculate the tree admittance | 

products and accumulate the values of the 
terms with same degree so that V(Y) or 


WW, rsa or & U (¥) will be obtained. 





designated as "W(L,J)." 
Ye aie inerereeal — —_ 


Je oe 


| 
| 
(3) Store these results in matrix form and | 
| 
| 





Block F ' 
— oe a wsstan apie Lo 


it Call. Mioratietne UResult": 
| 





Form the denominator and the numerator of the 
network function in polynomials. 

(2) Reduce the degrees of "S" if there are 
common factors in the numerator and 
denominator. 


2 eS) Print out the network eens) 


Block G 


| Prin Print out: program time | 
a ee 


ts more than one problem is 








inserted in the program. 
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4.6 Topological Formulas. 


The topological formulas [6]-[8] used in this program are listed 





below: 
Vi) 
Yap = Wye (4-6-1) 
VA i (4-6-2) 
EL Wy (1) + We eX) — Y)-—Wy a (1) 
Les lh a A AS on a (4-6-3) 
> VCY) 
Z.2= a (4-6-4) 
a ae (4-6-5) 
| 
; TO ag gee a came ZN 2 < pe (4-6-6) 
2. 
Nah age 
a (4-6-7) 
} 
ALW, 2(¥) + We, v0Y) — W201) —Wyy¥) 
Mh a = LW,z ) 2,1 U2 152 }) (4-6-8) 
ec Wh, CY) 
1 
x] Wa CY) + Wye YY) — Wi 2 OO) = Wier OY) 
A, = 22 00 + Mb 0) Swe O—vhie I] = =< J) (4-6-9) 
We, 2 t¥) 
SLM 2 CO Wy 270) — WA 2 OW (Y) 
A,,= BLS OO PO BM era ib _ = (4~6-10) 
We 2(Y) ‘: 
we W,' 5e = wi r( —We, (YT) 
fe EN OREO SO Mat be eins 


W,, 4 (yY) 
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5. Conclusion 

The results of the two new tree~finding programs show that the 
computation times and consequently the execution times for the deter- : 
mination of the network functions are greatly reduced compared to 
other methods. Furthermore, the program for determining network 
functions in Appendix LIT can be modified for generalized network 


analysis and synthesis, 
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APPENDIX I 
The digital computer program for finding trees by using the 
dichotomy method, as described in section 3.1, is presented below. 
Language: FORTRAN IV 
Machine: IBM 360/67 
Problems: 7 


The linear graphs corresponding to each problem are shown below. 









“a” & 





(f) Problem 6 


(g) Problem 7 


PRINT CANDITION="T2//, 
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et Wax & oe = web 
= a a wy ic CZ wt & o 
Qn & qf UU (I> «© Y Fr F& - 
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CO GT wet Azer & ke Ww Co foal 
TO rf > Tete ZO Ww FU OC = od 
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UOUO 


BLOCK 


9KSyKAgLTMP,CONTL »NX,MTREEF) 


0 TQ 300 
KEY 
GO 


+ 

0 

E 
(CONTL(M).EQ.0.0) GO TO 100 


COMUUSTRUCK HEL 
COZ ZmUOmMHOZS 





Y »KSpKA,LTMPypCONTL ¢NX_MTREF) 
yKAyLTMPyCONTL »NXyMTREF) 
KAgLTMPyCONTL »NXyMTPEE) 
KAyLTMPyCONTL 9 NX »MTREE) 

EY 9 KS,KAghL TMP, CONTL y NX» MTREE) 
KAyLTMPyCONTL y NX yMTREE) 


TIM MB KEV yg KS eKAeLTMPyCONTL pNXeMTREE) 
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APPENDIX II 
The digital computer program for finding trees by using the 
» T-triangle method, as described in section 3.3, is presented below. 
Language: FORTRAN IV 
Machine: IBM 360/67 
Problems: 7 
The linear graphs correspondin g to each problem are the same 


as shown in Appendix I. 
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NODES= 4, EDGES= 6, PRINT CONDITION= 1 
EDGES CALCULATED= 6 

INPUT LINK 

te 


INPUT T-TRIANGLE 


2 0 
3 4 6 
TREES 
a a 
1 2 4 
am & F 
eS  ¢€ 6 
i 3 6 
1 4 6 
L 3 6 
we of 4 
eos 3 
2 6 4 
z# 6 53 
3 4 6 
* 2 .% 
TOTAL NO. OF TREES= 13 


TIME= 0.63 SECONDS 


NODES= 5,» EDGES= 7, PRINT CONDITION= 1 
EDGES CALCULATED= 7 

INPUT LINK 

rok 2. o 


INPUT T-TRIANGLE 


TOT 
a 
~! 


WEOABIVIUUE ES RWRWWWWRRNN WW OW 
PAQW ANA ANNA OUI SUID 
NEHA NAWARWOP HANAN AO NANO I 


1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
2 
2 
2 
2 
2 
2 
2 
2 


TOTAL NO. OF TREES= re} 
TIME= 0.06 SECONDS 








NODES= 4, EDGES= 6, PRINT CONDITION= 0 
EOGES CALCULATED= 

INPUT LINK 

ti 


INPUT T-TRIANGLE 


3 5 6 
TREES HAVE BEEN FOUND, BUT NOT PRINTED 
TOTAL NO. OF TREES= 16 
TIME= 0.06 SECONDS 


zor 


NODES= 5, EDGES= 10, PRINT CONDITION= 0 
EOGES CALCULATED= 10 


INPUT LINK 

t Se a eh 

eS J 4 

x 

1 
INPUT T=TRIANGLE 

1 

2 5 

3 6 & 

4 7 910 

TREES HAVE BEEN FOUND, BUT NOT PRINTED 
TOTAL NO. OF TREES= 125 


TIME= 0.13 SECONDS 








NODES= 6, EDGES= 15, PRINT CONDITION= 0 
EDGES CALCULATED= 15 
INPUT LINK 


ee 


or 2 


Pelee 
p= 


= 


MPWNe 2 eter 
v0 
Cc 
mr 


T-TRIANGLE 


ad ©] 


1 

11 13 

12 14 15 

TREES HAVE BEEN FOUND, BUT NOT PRINTED 
TOTAL NO. OF TREES= 1296 


TIME= 0.38 SECONDS 


OOId 
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NODES= 7» EOGES= 21, PRINT CONDITION= 0 
EOGES CALCULATED= 21 


INPUT LINK 

i ae fs OF 
» te iE Ff 

s+ L ff fj 

EL f 

kt 

1 
INPUT T-TRIANGLE 
1 

fy i 

a 8 32 

& 39 13 N6 

5 10 14 17 19 

6 11 15 18 20 21 


SOT 


TREES HAVE BEEN FOUND, BUT NOT PRINTED 
TOTAL NO. OF TREES= 16807 
TIME= 3.00 SECONDS 
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NODES= 8, EDGES= 28, PRINT CONDITION= 0 
EDGES CALCULATED= 28 


INPUT LINK 

teem Le tebed 
pono th b sb 

i Si Se ae | 

» © 2 

te t 

Lt 

1 
INPUT T-TRIANGLE 

1 

2 8 

3.9 14 

4 10 15 19 

5 11 16 20 23 

6 12 17 21 24 26 
713 18 22 25 27 28 


TREES HAVE BEEN FOUND, BUT NOT PRINTED 
TOTAL NO. OF TREES= 262144 
TIME= 32.87 SECONDS 


je 


APPENDIX ILI 
A computer program for the determination of network functions. 
The same networks as given in examples2.1 and 2.3 are solved in this 
program. 
Language: FORTRAN IV 
Machine: IBM 360/67 


Problems: 2 
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PRINT 29 
29 FORMAT €//15X, "ORIGINAL T-TRIANGLE*/) 
DO 26 I=1,NV 
DO 26 J=1sI 
26 NT(1,TpJI=NK(T oJ) 
00 36 f° =1,NV 
30 PRINT 31, (NTCL,I,d) pJ=lel) 
2 FORMAT (/15X,2513) 
C BLOCK C 
N3=3 
NG=4 
IF (IGR«EQ.1) N4=3 
ME=B 
: IF (IPORT.EQ.1) MF=2 
C BLOCK D 
CO $9 $1470 13914915916917 918 Io 
? ? g 9 % ? ? g 
ll NI=N 
GO TO 100 
CocccccsecccccceaesTO YIELD THE NET FOR Wel! 
12° [K2N3-1 
JK=1 
60 TO 200 
Gisvevews! SeeaeccevesseulG YIELD THE NET FOR W2,2! 
13° 1K=N4-1 
JK=2 
60 To 200 
Cutscene weseoossenld YIELD THE NET FOR W1,2 
14° tksi 
JK=1 
60 To 200 
Gowns fRisavevesevelO YIELD THE NET FOR wWi',2¢ 
if IK=N4-1 
JK=3 
IF(IK.LT.JK) GO TO 41 
GO TO 200 
G1 PRINT 445L,1K, UK 
44 FORMAT (7/720Xo" L=",1295X_"IK="gI12,y! 
DO 42 I=1eITERM 
42 WV(L,I)=0. 
PRINT 43,0WV(L oJ), J=Ly I TERM) 





JK=",12//) 


WV(LeJd=",15F6.2) 


(/10X,*NO TREE*//* 
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l 
esccoeeetO YIELD THE NET FOR W 


1,2! 


eoeeeoeTO YIELD THE NET FOR W2,1! 


AND1,1! 
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MODNET( IK, JKaNoh y NK NT NI yg TPP ye Ve ISy TPPNyLABEL gL AL yN4) 


ovsvecealS YIELD THE NET FOR U--SHORT 2 
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VLIV2 Sm rn rr earn 
W2,2°(Y) 
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PRINT OUT THE NETWORK FUNCTIONS 
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1.00S** 0+ 4.00S** 1+ 10.00SKe 24 12.00S** 34 6.00S** 4 
1.00S** 0+ 4,00S** 1+ 10.00S#« 24 12.00S** 3+ 6.00S** 4 
1.00S** Ot 7.00S** 14 12.00SKe 24 6.00S** 3 


TIME= 0.25 SECONDS 











INPUT DATA 


EDGES= 5 GR CONDITION= 2, 2-PORT 


NODE 


IPP(K) 


ted ted 


SNA 


Nats 


NX NY NK(NY-12NX)=K 


ol et OA) 


2) 
UJ 
Qa 
oO 
z 
z 
S 
a 
oO 
z 
~ 
4 
2 
z 
—_ 
= 
uw 
oO 
2 
| 
~~ 
Oo 
z 
Ww 
Q 
oO 
= 
~ 
oO 
z 
Ww 
o 
(=) 
{re 
~ 


ORIGINAL TF-TRIANGLE 


v(y) 


MINED IAAT 
NNN SIN 


et etedt ed HON 


WWW WW WL 
UWI WU LL) 
Ceerrrre 
a 


OnOANHON 
! 
He He HH He eH 


oooo0ocsdo 
eevee oes 
ANAN St tNN 


ee 
I, eR, AI, AI Ee, A, I 
>>> >> >>> 


Qooo0ceca 


400400008 


UUW a uu Wu 
WWW 
Cervera 
aa 


OF TREES= 
0.0 


TOTAL NQ. 


6.00 4.00 1.00 0.90 


2-00 


0.0 


WV(L eID= 


Wl,l'(Y) 





0.9 
0.0 


1.00 


2.00 4.00 
1.00 0.0 


3.00 


4.00 
2-00 


0.0 


OF TREES= 
0.0 


OF TREES 


0.0 
0.0 


TOTAL NQ. 


TOTAL NO. 
WV{L I) 


W(L eI) 


W292"(Y) 
W1 92Y¥) 





0.0 


[2] 
e 
i=) 
MS Inn sf Ateneet 
o 
NNN OS e MMOMARS 
=~ 
mn 
LULU LU LU) o LRU UU 
WW i} i] WU LU 
ae oa dad * ww [= qo Gi. dos 4 
EREEE na 9 ray ay Say ey 
(2) 
o ™ 
" e 
RMOHOnm HW + I atHAON 
1 | Ww ~~ 1 
HHEHH WwW el RRR HH 
RRB H cx HHH 
NNMMY —- © NAMM 
©oooog e Ooo0eo 
@Go0o0o0 th O Oo00o 
Q0000 0 oo000 
ooocoo un oo°co0o 
eeeee e e@eerve 
NNN O OC I Netrnedt 
= © cl 
oO 
wd 
ec 
HR wet ke U wnnnn 
mamamenm OC = a ey te 
>>>>> Fe >>>>> 
mee St Cae? et - ~_ er te ee eee 
O0000 1 ~- O0A00 
Qo00c0 ~ - oOoceoc 
Cxxrrae > e taeter 
40a.0,.0.0. » 4 N Oo28002080. 
- 
WLU LL ue = Wu 
LY OL LL) Lk ea WIL 
exeaorer =z Cerrarer 


EReEER aa 


131 


OF TREES= 
0.0 


TOTAL NO. 
0.0 


WV(L,LT)= 


2.00 2.00 1.00 0.0 


2.00 


-™ 
> 
— 
oo 
N 
e 
~ 
=z 


om eee 


OF TREES= 
0.0 


TOTAL NO. 
0.0 


WV(L I )= 


4.00 2-00 0.0 0.0 


2-00 


W2,1°CY) 


5 
1.00 


OF TREES 
0.0 


TOTAL NO. 


WV(L »DD= 


1.00 0.0 


0-00 3.00 


0.0 


UC Y) 


ANON 


TOTAL NO. OF TREES= 
0.0 


WV(L al )= 


0.0 


2.00 2.00 1.00 0.0 


0.0 





eel 


Z12=221 


Z22 


Y1l 


Yi2=Y2l 


Y¥22 


V2/vl 


Vi/V2 


I2/11 


2.00S#* 
4.00S#* 
4.00S#* 
2.00S** 
-2.00S** 
2.00S*# 
2.00S** 
2.00S** 
2.00S** 
2.00S#* 
2.00S*# 
2.00S#* 
4. 00S#* 
2.00S** 
-2.00S** 
4.00S#* 
-2.00S#* 
2.00S#* 
2.00S** 
2.00S** 


PRINT 

Ot 6.00S** 
Ot 2.00S** 
O+t 2.00S** 
O+ 6.00S** 
0+ 1.00S** 
Ot 6.00S** 
Ot 3.00S** 
Ot 6.00S** 
O+ 3.00S** 
Ot 2.00S** 
O -1.00S** 
Ot 2.00S** 
Ot 2.-00S** 
Ot 2.00S** 
O+ 1-00S** 
0+ 2.00S** 
Ot 1.00S** 
O+ 3.00S** 
0 -1.00S** 
Ot 3.00S** 


4.00S** 
4.00S** 
4.00S** 


4.00S** 
1.00S** 
4.008** 


1.00S** 2 


1.00S** 
4.90S** 


1.00S** 


2 


OUT THE NETWORK FUNCTIONS 


1.00S** 
1.00S** 
1.00S** 


Oe ee ee oe ee ee ee a ee ee ww we ew eee we ee wee em em ee we eee eee ee ee ee owe 
ee ee a ee eee a we ee ee ew ee we © ae ee wm mem ew ew we we eee ew ww wwe wm = 
ee oe ne ee ee ee ee we ee ee ae ee wee ew ee ee ew me ew we we ee eee ee 
Oe ee ae ee ae ee ae ee ee ee ee ee ee ee ee ee re ee a em ee OM em ee ee we eo ee ne ewe eee 
ee ee ee ee ee ee ee ee ee re ew we Ow ee ee ew eee Be ew ww ew we een 
ew we we we a we Oe wm oe ww we wr we ow we Oe ee ew wee SO ee 
Se et ee mes Se ee ee OS Oe Oe ee Se Sew eee ee were eee ee eee 
Pe mee 6 em ee wwe we em Oe we ee re ee ew ww we Oe ee ow ee ew we ee Bw ewe ew we wee oe 
2 3Pe eee ee we ee Oem OO OO Oo wr 8 ee 2 OO Oe ewe Be ee ew wee wee ee ee = 


See BOO OS MK ee ee OOOO eS ee eS ew we eM Se ee we ee SS ee Oe Bee Swe ee ew wwe Se wee eee ee 


a = ee 


*» 


SGNOIJ3S O8°O 


€ #4S00°T +2 **SO00°Y +T **SO0°C +0 *#S00°Y 


es ee ce ee a Ow ae eee a ee ee oe ee ee we eS ee ee eee ween ee = 


cf **SOO°T— GO **S00°2 





=3wIl 


= 21/1! 
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